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Abstract 

The effect of quenched random ferromagnetic bonds on the antiferromagnetic 
correlation length, of a two-dimensional Heisneberg model is studied, 
applying the renormalization group method to the classical non-linear sigma 
model with quenched random dipole moments. It is found that the antiferro- 
magnetic long range order is destroyed for any non-zero concentration, x, of 
the dipolar defects, even at zero temperature. Below a line T oc x, where T is 
the temperature, is independent of T, and decreases exponentially with 
x. At higher temperatures, it decays exponentially with yo| ff /T, with an effec- 
tive stiffness constant which decreases with x/T. The results are used to 
estimate the three-dimensional Neel temperature, which decays linearly with 
x at small concentrations, and drops precipitously at a critical concentration. 
These predictions are compared with experiments on doped copper oxides, 
and are shown to reproduce successfully some of the prominent features of 
the data. 
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I. INTRODUCTION 

Consider a Heisenberg antiferromagnet, with a concentration x of quenched random near- 
est neighbor ferromagnetic (FM) bonds. These bonds compete with the antiferromagnetic 
(AFM) order, and introduce frustration into the problem. If the FM exchange on the "im- 
purity" bonds is sufficiently strong, then the two spins at the end of each such bond prefer 
energetically to be parallel to each other, and perpendicular to the background AFM stag- 
gered moment. The staggered moments on the other sites then cant in this perpendicular 
direction, and at large distance this canting angle decays with distance r as l/r d_1 in d di- 
mensions, similarly to the decay of magnetic moments in the presence of a magnetic dipole. 
|Ij] This follows from a mapping of the low temperature equations for the spin configuration 
at the minimal energy onto the Laplace equation. At sufficiently low d such impurities de- 
stroy the AFM long range order at all finite temperatures T, giving rise to a spin glass like 
phase. The present paper presents a renormalization group (RG) analysis of the AFM 
correlations in the presence of such quenched randomness. As we show, this complex ran- 
dom problem is exactly renormalizable in two dimensions (2D), allowing a detailed study 
of the dependence of the 2D AFM correlation length £ 2 d on T and on x. This also allows 
an e-expansion in d = 2 + e dimensions. The three dimensional Neel temperature T^(x) of 
lamellar system is estimated by a model of weakly coupled planes. 

It is hoped that our model can describe the T — x phase diagram of various antifer- 
romagnets with random FM bonds. One motivation of the present study arises from its 
possible relevance to the understanding of the doping dependence of the magnetic order in 
the lamellar copper oxides. Experimentally, doping such oxides leads to a rapid decrease in 
both £ 2 d and Tjy. PHlQj Experiments on doped La2Cu04 show that above x = x c ~ 2% £ 2 d 
remains finite even at zero temperature, and there is no AFM order. This strong effect of 
the doping has been attributed to frustration, due to strong FM exchange on the Cu-O-Cu 
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bonds which have localized holes due to the doping. This frustration was also predicted 
to yield a magnetic spin glass phase for x > x c , as recently confirmed in detail in doped 
La2Cu04. [Tl[] The experimental verification of this spin freezing above x c [|i||S||T0|-|r3| con- 



firms the picture of localized holes. This localization is also confirmed by direct conductance 
measurements at x < 5% and low temperature T. P|.|TZ[ 

The description of doped lamellar cuprates by quenched FM bonds requires several as- 
sumptions, which will be discussed in detail below. In particular, at high temperature T the 
"dipolar" moments which describe the dopant bonds may fluctuate, turning this aspect of 
the problem into one which requires a combined annealed and quenched averaging. Indeed, 
Glazman and Ioselevich (GI) studied this problem in its annealed limit, and calculated 
£ 2 d to leading order in x/T (See also Ref. |15|].) Although the locations of the dipole-like 



impurities are randomly quenched, each impurity involves an effective dipole moment m(r^) 
which is still free to reach annealed equilibrium in the presence of all the other dipoles. At 
low temperature, these moments freeze as in a dipole glass. However, unlike the dipole 
glass, where the interactions are fixed, the interactions among the dipoles are mediated by 
the AFM spin background, whose behavior also depends on temperature, concentration and 
configuration of the dipole moments. In the absence of a simple systematic way to handle 
such a combined quenched-annealed problem, GI stopped their explicit calculations at the 
low-x expansion. Here we argue that at sufficiently low T the dipole moments, which inter- 
act via randomly quenched dipole-dipole interactions, freeze in a random spin glassy way. 
T6| , |T7H At low temperatures we thus compare the experiments with our quenched theory. 



The actual fitting of data from the cuprates should involve some interpolation between the 
annealed and quenched limits. As we show, the annealed and quenched calculations coincide 
to leading order in x/T, and therefore our quenched results supply a good interpolation over 
the whole range. Other assumptions, concerning e. g. the mobility of the holes forming the 
dopant bonds, will be discussed below. In any case, our model should give a good descrip- 
tion of the lamellar cuprates at low T, and should describe other lamellar antiferromagnets 
doped with FM bonds. 



A second major motivation concerns the fact that, as we show below, the quenched 
random dipoles are coupled to the gradient of the order parameter, and therefore they are 
equivalent to correlated random fields, whose correlations in momentum space are propor- 
tional to the square of the momentum. Such fields lower the critical dimensions of the 
random field TV-component spin model by 2, from 6 to 4 for the upper critical dimension, 
and from 4 to 2 for the lower one. In Ref. |18| this has been established for the limit Af — > oo. 



Here we show that the lower critical dimension is shifted down from d = 4 to d = 2, for all 
M > 2. As a result, both the temperature and the variance of the random dipole moments 
(which is proportional to the concentration, x) are marginal (in the RG sense) at d — 2, 
allowing for an analytical solution of the recursion relations. This marginality of the ran- 
domness is related to the 2D infrared divergence of the Villain canted states. |IJ Within the 
one-loop approximation, we obtain an exact expression for the exponential part of the 2D 
correlation length, which remains finite at all non-zero x. 

We describe the system by the reduced Hamiltonian (i.e., the Hamiltonian divided by 
the temperature T) 

7~t = TYpure + ^intj (1) 

where 7i pure is the classical non-linear sigma model (NLcrM) for the pure (non-random) 
system, representing the long wave length Hamiltonian related to the fluctuations of the 
unit vector n(r) of antiferromagnetism, 

ftpure = 7T t J dr E(^%) 2 > * = T /Pf ( 2 ) 

Here i = 1, .., d and \x = 1, .., Af run over the spatial Cartesian components and over the spin 
components, respectively, p s is the stiffness constant and <9j = d/dxi. This Hamiltonian is 
known to give an excellent description of the undoped antiferromagnet, both theoretically 
fl9|| and experimentally |4j]. The success of this description results from the fact that, 
although the problem involves quantum spin fluctuations, these can be integrated out at 
any finite T, causing only a renormalization of p s . 



Hint is constructed [Tj] to reproduce the dipolar effects at long distances: Denoting by 
a(r^) the unit vector directed along the doped bond at rg, and by Mm(r£) the corresponding 
dipole moment (where m(r^) is a unit vector giving the direction of the dipole, and M is its 
magnitude), 

Wfart = T / rfr H f i( r ) ' <^ n ' ( 3 ) 

t J i 

with 

fi(r) = M]T<5(r-r,)a,(r £ )m(r,). (4) 

Note that since n is a unit vector, diXi is perpendicular to n and hence 7ii nt contains only 
the M — 1 components of the TV-component vector fj which are transverse to n. However, 
since the vector n varies with r, all the components of m may enter at the end. As stated, 
GI treated the variables m as annealed. Here we treat all the variables r^, a(r^), and m(r^) 
as quenched. Denoting quenched averages by [...], we write 

[a(r)] = 0, 

K(r)] = 0, 
[ai(r)aj(r')] = <%<!>(r - r')x/d, 
[m^{r)m v {r')\ = Q8^8{r - r'), (5) 



so that [fifj] = and 



[U(r)f Jl/ (r')]=X5, I/ 5 lJ 5(r-v'), (6) 



with 



A = M 2 Qx/d = Ax, (7) 

where A = 0(1). At low T one expects Q ~ 

One can see now that 7i m t represents random fields with quenched correlations: Fourier 
transforming the variables in Eq. (||), 7imt can be written in the form 



Hint = Wfi I ^EM^KC-k), (8) 

with the random field h(k), 

h^) = iYk 3 jdTf m {r)^\ (9) 

3 

which has quenched correlations 

[h,(k)hl(k f )} = \k\ u 5(k-k'). (10) 

Such correlations shift the critical dimension of the random field Heisenberg problem down 
to 2. A heuristic way to show this follows Imry and Ma [[2(J in assuming an ordered state 
and considering the transverse spin fluctuations, JA ± (t). In momentum space, ■M J ~(q) = 
G ± (q)h ± (q), with G ± (q) ~ 1/q 2 , and /^(q) denoting the transverse random field. Thus 

[M ± (r)M ± (r')] = (^-) 2d ( d d qd d q' x 

G i (q)G ± (q')[/ i ± (q)^(q / )]e i(q - r+q '- r ' ) . (11) 

For [/i /i (q)/i I /(q')] oc 5 f±u q e 5(q + q'), this integral diverges for d < 4 — 9, implying that the 
assumption of long range order is invalid. This identifies the lower critical dimension as 
di = 4 — 9, and for our case, d% = 2. The calculations below support this picture. In 
addition to giving an exact solution for the problem at hand, we note that the present 
formalism might also be used as a starting point for a double expansion, in e and in 2 — 9, 
aiming at other random field problems. 

The conventional method to treat the NLcxM expands the order parameter unit vector 
n about a spatially uniform ordered state. pT] , p2| Using replicas to handle the quenched 
randomness we have found that this approach generated local random fields, which seem to 
break the symmetries of the original model. Similar problems were found for other problems 
near 4D, where they required a resummation of the perturbation expansion. [F23J In our case, 



the problems may have come either from the replicas or from the assumption of a spatially 
uniform state. Such a state does not exist when there is some local freezing of the moments 
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in random directions (as happens in the random case discussed here). We avoid both of 
these by going back to the original RG approach by Polyakov, and by treating the 
randomness without replicas. In 2D, this allows us to obtain for all values of p s x/T. 

We show that the quenched dipoles suppress the antiferromagnetic correlations, so that 
the correlation length is a decreasing function of p s x/T, remaining finite at any non-zero 
x, even at T = 0. This implies the destruction of the antiferromagnetic long range order in 
2D for any concentration. As x increases, ^2d{T = 0) decays exponentially, representing the 
sizes of the Imry-Ma domains for this system. 

In order to compare the results of our model with experiments performed on doped 
lamellar cuprates, we consider the 3D ordering of a system of weakly coupled planes by 
using the relation a£% D ~ 1. Here, a represents either the interplane coupling, that is the 
relative interplane exchange, J±/J ~ J±/27cp s , or the in-plane relative spin anisotropy (in 
the presence of which even an infinitesimal coupling suffices to yield 3D ordering). Although 
it is not expected to give the correct 3D critical behavior, this procedure proved to give 
excellent results for the 3D Neel temperature in the pure case. [|J This "mean field" procedure 
is also justified by an RG argument: to linear order in a, the RG recursion relation for 
a is a' = e 2£ a, where e e is the length rescale factor. After t iterations of the RG the 
effective coupling between planes involves renormalized spins, contained in the area e 21 of 
the renormalized cell. A 3D behavior is expected when this effective coupling becomes 
comparable to 1. As we show in Appendix B, similar results for the phase diagram are also 
found from integrating the RG recursion relations in d = 2 + e dimensions. We obtain an 
explicit form for T^(x), and compare it in detail with available data. Our results reproduce 
prominent features of the observed phase diagram, in particular the fast decrease of the Neel 
temperature with increasing x and the disappearance of the 3D long range order at x ~ 2%. 

The outline of the paper is as follows. Section II discusses the RG procedure, and Sec. Ill 
describes the RG recursion relations for the quenched averaging. The 2D recursion relations 
are then solved in Sec. IV, and the resulting £ 2 d is used for estimating the 3D phase 
transition line Tn(x) in Sec. V. Section VI then contains a discussion of the alternative 
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annealed averaging. The results are compared with experiments on doped cuprates in Sec. 
VII, and discussed in Sec. VIII. Details of the calculations and extensions to d = 2 + e are 
given in the Appendices. 



II. THE RENORMALIZATION GROUP PROCEDURE 



Following the RG approach of Polyakov, P^ , p5| we decompose n(r) into a slowly varying 



part, given by the unit vector n(r), and M — 1 fast variables At (r), such that 



n(r) = n(r)yi - <f) 2 (r) + ^ M (r)e M (r), 

2 W= (12) 

The unit vectors ii(r) and e M (r), /i = l,..,J\f — 1, form an orthonormal basis. The Fourier 
transform of the fast variables 

M (r) = (2vr)^|^ r M (q), (13) 

is restricted to wave vectors q in the range 6" 1 < q < 1. The upper bound is the inverse 
of the microscopic length (which is measured in units of the lattice constant), and b is the 
length rescale factor for the renormalization procedure. These q values are to be integrated 
out. After the iteration the correlation length £ is renormalized into £/b. 

The Hamiltonian 7i requires the derivatives of n(r). Using the relations n • <9jn = 0, 
n • e M = 0, and e M • e u = 5^ u , we set 

= ^ £f e M , d ie , = £ Afe, - 5f n, (14) 

where A 1 ^ = e M • die u = —Af w . Then we find 

JV-l 



d i n = n{d l ^/l-^- ]T 

n=i 

+ X 1 e At {^ + J B,r v / i r ^+ EA^}- ( 15 ) 



s 



We show in Appendix A that the functions A^ u , which give the first derivatives of the base 
vectors e M , can be eliminated by a suitable gauge transformation when one ignores higher 



order derivatives. 24,2c Therefore, these are omitted in the following. 



In terms of the new variables, the Hamiltonian TC, to order 2 , reads 

n = H +H 1 +H 2 +H 3 + H i , (16) 



with 



and 



H ° = YtJ dr T,{mr)f + (^0,(r)) 2 }, (17) 

ifi 

n i = Yt[ *E5f(r)Br(r){^(r)^(r) - <W 2 (r)}, (19) 

iufi 

n s = -- j dr]T5f(r){^(r)^(r) + ^f(r)0 2 }, (20) 

ifi 



^ = 7 / *EK( r )^M(f) - ^(rWj(r)}. (21) 
Here we have introduced the notations 

« 4 (r) = n(r).^(r), ^f(r) = e,(r) ■ f,(r), (22) 

for the longitudinal and transverse components, respectively, of f, in the new variables. 
Indeed, Eq. flT8|) represents the bare form of TC- m t in this system. In TCo the slow variables, 
.Bf , are separated from the fast ones [One should notice that (d^) 2 also yields the 
contribution (1/t) J dr B% di<p M . However, this term vanishes upon Fourier transforming, 
as </> M pertains to the large q portion of the Brillouin zone, while the slow variables belong 
to the small q values, q < b^ 1 .} The Hamiltonians 7^2, 7"^3, and TC^, of order 0(B 2 ), 0{B 1 ) 1 
and O(B ), respectively, will be treated in perturbation theory. 



III. RENORMALIZATION GROUP EQUATIONS 



The first step in deriving the recursion relations involves integration over the fast variables 
M . This requires the Green's functions, 

(^(r)0 l/ (r')) = ^G(r-r') ) 

G(r) = (2vr)- d / dqe^G(q), G(q)=t/q 2 , (23) 

where (..) denotes a thermal average with 7io- As we shall see below, to leading order in 
e = d — 2 we need G(r) only at strictly 2D, where 

G(0) = i-ln&. (24) 

Hence (0 2 ) is small for In 6 -C 27r/t. In practice, G(r) is significantly different from zero only 
for 1 < r < b, where it is approximately given by the 2D Coulomb interaction 

G(r)«^m-. (25) 
2tt r 

We next turn to the perturbation expansion in TL 2 , H3, and TC±. The first order yields 
= (H 2 + H 3 ) = 1 / dr£(Ef (r)) 2 (2 - AT)G(O) 
-^/dr^Sf(r)^(r)(jV-l)G(0). (26) 



t/Li 



The first term here represents the leading order renormalization of 1/t, as usual. [P2,P1| The 
second term, which is linear in Bf, is similar to the initial 7i- m t, or to the equivalent Eq. 
(Pf). In fact, this term contributes to the renormalization of the transverse components of 
f,. 

Higher order perturbations contain higher powers of 0, which yield higher powers of G 
and hence of t. In the following, we keep only leading powers of t. Neglecting the terms 
involving products of two G's, the second order perturbation yields 
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+ i| dndr 2 E^( r i) u *( r i)^( r 2) 8 ^ G! ( r i2) 



'J A' 

' f rfr lC /r 2 ]T J Bf(r 1 ) J Bj t (r 2 )n l (r 1 )^(r 2 )G'(r 12 ) 



2t 2 



r i2 = ri - r 2 . (27) 

The first term here is -B-independent. In principle, it gives rise to an interaction between 
the dipoles: Inserting G(r) and the explicit expressions for gff [Eqs. (^) and (E^)] we rewrite 
this term in the form 

ftdd = - hem±(r k ) ■ m±(r i ), (28) 
t u 

where mi (r^) is the component of the dipole moment at which is perpendicular to n, and 

1 „ 1 

f o ( a K0 -r H )(a(r £ ) -r ke ) 

{2 2 a(r fe ) • a(r £ )|, (29) 

with < 6. Apart from trivial factors, this reproduces the effective dipole-dipole interac- 
tion found in Ref. |14| |. There, Eq. (|28| ) was used to integrate over the variables mi, treating 



them as annealed variables. In the present calculation we treat the dipoles as quenched, and 
therefore Eq. ( p!jj) simply represents an additional constant to the energy. We return to this 
point in the following. The other two (5-dependent) terms in the second order perturbation 
Hamiltonian, Eq. fl2"7p, will contribute to the renormalization of the temperature and the 
variance of the dipolar quenched interaction. 

Finally, the third order perturbation Hamiltonian, keeping terms up to order B 2 , is 

H {3) = \{Hf) + l((HlH 3 ) + (H4HI) + (H 2 Hl)). (30) 



6 N *' 2 

Integrating out the variables <ft, it is seen that the first term here is independent of B\ 



it contributes further to the dipole-dipole interaction. The next term in fl3"U| ) yields an 
expression linear in B, 
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" ijk A* 

{^^(r 1 )^(r 2 )^(r 3 )9 2j 93 fc G(r 21 )G(r 31 ) 

V 

-2« i (r 1 )^(r 2 )« fc (r 3 )9 2 ^3 fe G(r 23 )G'(r3i)}, (31) 
which again contributes to 7"^nt, while the last term there gives a B 2 contribution 
l -(H A Hl + H 2 Hl) = ^J dv x dv 2 dv^Y. 



x 

ijk M 



{- J B < "(r 1 )S i ''(r 2 ) Ui (r 1 )u i (r 2 )u fc (r3)^ji,Gf(r 13 )G?(r 2 3) 
+2^5f(r 1 ) J BJ(r 2 ) Mj (r 1 )^(r 2 )^(r 3 )9 3fe G(r 12 )G(r 23 ) 

+ E(^(ri)fir(ri)^(r 2 )^(r 3 )- 

V 

5r(r 1 ) J Bf(r 1 )^(r 2 )^(r 3 ))a 2 ^ 3fc G(r 12 )G(r 13 )}. (32) 



Except for the first term in Eq. (p6|), all the generated terms involve the longitudinal and 
transverse components of the vecotrs f,, and , Eq. (^2j), which depend on the quenched 
random variables r^, a(r^), and m(r^). Using Eq. (|g) we thus find 

[#f ( r )#j ( r ')] = \8ijS^8(T - r), 
[«i( r )«j( r ')l = X5ijS(r-r'), 

K(r)^(r')] = 0. (33) 

We now obtain the recursion relations of the RG. Consider first the quenched averages of 
the integrated Hamiltonians Tt^\£ = 1, 2, 3. This will constitute the renormalization of 1/t. 
Rescaling the lengths by b^ 1 and the slow derivatives B% (q) by 6 , the new temperature 
prefactor, multiplying the integral over (-Bf ) 2 , obeys the RG equation 



- = b d - 2 
f 



1 2-N, . 1 - Af \ 



In 6 + Info, (34) 

t 2ir 2ir t J' y J 

which is valid to first order in e = d — 2, t and A. To obtain this equation, we have used 

Eqs. ( |33"1) and the relation 

£ / dr 2 (d 2j G(r 12 )) 2 = tG(0). (35) 

3 
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The terms in the Hamiltonian linear in B- 1 remain as quenched random contributions. 
They renormalize /j M and yield a renormalization of its variance A. To obtain the recursion 
relation for A, we collect all terms linear in B and write them in the form 



(36) 



with 



Tffr 



flf(r)(l-^(jV-l)G(0)) 

+ 7 / dr 1 ^2u i (r)g l j(r 1 )d lj G(r - ri) 
t J . 

jk 

{Y^g^^g^g^d^Giv, - r)G(r 2 - r) 

V 

-2ui(r)g^(r 1 )u k (r 2 )d lj d 2 kG(r 1 - r 2 )G(r 2 - r)}. 



(37) 



We then find the variance of Tf(r) by the one-loop calculation, which to leading order, using 
Eqs. (pi), yields 



rf(r)rf(r')] =5 ll ,5^5(v-v')x 
\(l-{N- 1)G(0) + -(2 - AT)G(O)) 



(38) 



Hence, the recursion relation for A is 



b 7 



t(Af - 1) + A (A/" - 2) 



2tt 



Info 



(39) 



One also needs to consider the fluctuations of the random terms around their quenched 
averages, as well as new terms, which were not included in the initial TC, but are generated 
by the renormalization procedure. However, these are irrelevant. Let us take as an example 
the last term in which we may write as 

1 



2t 



dv l dv 2 Y J W l3 {v l v 2 )B^{v l )B>;{Y 2 ). 



(40) 



i-iij 
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Physically, this term describes a random interaction among the gradients of the order pa- 
rameter n, which are absent in the original problem. As the ensemble average contribution 
of this interaction has been analyzed above, we need to consider here the deviation 



5W l3 ( ri r 2 ) = 

-- Uitr^Ujfa) - X5 ij 5{r 1 - r 2 ) G(r X2 ) 
Using Eqs. (H), it is easy to see that 



X 2 

— [SwSjySfa - t[)5(t2 - r' 2 ) + 

V 

5 if 5 ji/ 5(r 1 - r' 2 )5(r 2 - r / 1 )]G(ri 2 ) 



(41) 



(42) 



Since the range of G(r) is of order b, the correlations among these generated W*s are short 
range, and in practice the W's can be treated as uncorrelated, i.e. 

[5W ij {r l v 2 )5W iir {v' 1 r' 2 )\ = 
A[5 ii /5 7 y5(ri - r / 1 )<5(r 2 -t' 2 ) + 

SifSjffSin - r' 2 )5(r 2 - r[)]5(r 12 ). (43) 

A simple power counting then shows that (W/t) scales as b 2d ~ 2 , W scales as h d and hence 

dA 



-dA + 0(t,X). 



(44) 



Therefore, this generated random coupling is irrelevant. Similar arguments apply for the 
variances of all the other random terms which are generated in the renormalization proce- 
dure. 

We now follow standard procedures, and use an infinitesimal length rescale factor b = e Si . 



To linear order in e = d — 2, t and A, Eqs. ( p4|) and (|39|) now become 



d 1 _ 1 2-A/" 1 - jVA 
di~t~ e ~t + 2tt + 2tt T 
d L \_ \ l-N X 2-AT X 2 
di¥ ~ t t 2 ~ + 2?r t + 2vr 



(45) 
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dt 




M -2 2 


M-l 


~d£ ~ 


-et-t 




2?r 


dX 


-eX- 


A/"-3 x 










~di ~ 




f 2. 


2tt 



Combining these two equations yields 

\T - 2 „ A/" — 1 

(46) 

As noted above, d — 2 is the lower critical dimension for the NLcxM with quenched dipoles. 
Hence, the 2D problem is exactly renormalizable (as done in the next section), and one can 
obtain an e-expansion in d = 2 + e dimensions (as done in Appendix B). 

IV. THE CORRELATION LENGTH OF A 2D HEISENBERG SYSTEM 

We now proceed to calculate To this end we solve Eqs. (|46|) with the initial values 
t(£ ) = t , and X(£o) = A . The parameter £ represents some prefacing iterations. In the 
simplest case, we assume that £$ = 0, and thus that e io = L Q — 1. Other choices for L will 
be discussed below. The solution is particularly simple for the Heisenberg system, M = 3. 
At 2D one finds 

X(£) = X (l-^(£-£ ))-\ 

Both t(£) and X(£) flow away from the fixed point t = X = as £ increases. ^From the 
second of Eqs. ( f47f ) it is seen that Ao > t implies A > t, and vice versa. 
The standard scaling relation for the correlation length is 

£{t,\) = Jt(t(l),\(£)). (48) 

The correlation length £ is obtained from the matching condition at £ = £*, 

max(t(f ),A(f )) = 2tt, (49) 

where £* is chosen so that £(t(£*), X(£*)) is of the order of the renormalized lattice constant. 
In practice, this implies that at £ = £*,£ is a slowly varying function of its variables, which 
we denote by C(t, A). The first of Eqs. ( ^7|) gives 
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where X(£*) is equal to 2n for Ao > to, and is given by the solution of 

/ Ao \ 2 / 3 r/ Ao y 



Ao \ 2 / 3 |Y Ao \V 3 , Ao 



Ao 
2tt ! 



(51) 



for t > Xq. Equations (|48| ) and fl5"0D give the 2D correlation length, 



A) = L C(t, A)exp 



2tt 



1 



Ao 



L3A ^ A(£*) 

In the low temperature limit, Ao/to ^ 1) we stop iterating when A(f* 
consequently 



£ 2 d ~ ^oC(t, A) exp 



r 2vr 

3^ 



i 



Ao 
2tt 



(52) 
27r, and 

(53) 



It follows that ^2D is finite at any finite Ao, even as to approaches zero. This implies that 
at zero temperature the long-range order in 2D Heisenberg magnets is destroyed at any 
small amount of defects (as indeed predicted already by Villain. 0) This conclusion is also 
supported by the observation mentioned above, that Hmt represents correlated random fields. 
The exponential form of Eq. (|53|) is similar to that found for other random field problems 



at the lower critical dimension. |2(J Monte Carlo simulations §27j also suggest that the zero 
temperature correlation length is finite in 2D classical Heisenberg magnets with frustrated 
bonds. When the correlation length of the 2D system remains finite at zero temperature, it 
measures the size of the Imry-Ma domains, which is given by an exponential form like Eq. 



When Ao < t 1, we can approximate Ao/A(£*) by (1 — Ao/to) 3 [see Eq. (fTf]) ] and 
obtain 



/27T 

£ 2 d = L C(t,X)exp[ — 



1 Aq | Aq 
to 3to 



(54) 



The exponential part may be interpreted as a renormalization of the effective stiffness con- 
stant in the usual expression for the 2D Heisenberg model, 
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pf = p s (l-y + y 2 /3), 



y = A /to = AxPs/T, 



(55) 



where we have used Eq. (|7|). To leading order in xp s /T, this coincides with the expression 



which was obtained in Ref. |14| for an annealed system of dipoles. Indeed, up to the lowest 
order in X/t there is no difference between quenched and annealed averaging. 

Finally, we discuss the pre-exponential factor in the expressions for the corrletaion length. 
For A <C to, the prefactor C(t, A) ~ C(t ,0) is known: The two-loop [|19j and three-loop 
calculations, based on the quantum NLcrM, show that 

to 



C(t ,0) 



e_c_, 
8 2np s V 4tt 



(56) 



where c is the spin- wave velocity. 

At low tempreatures, we need the concentration dependence of the pre-exponential factor. 
This results from higher order loops: At 2D, t = and M = 3, the generalized recursion 
relation for A has the generic form 



^ = (3 2 \ 2 + /3 3 A 3 , 



with p 2 = 3/(27r) and with (3 3 of order /3|. The solution for this equation reads 



(57) 



( 



x 



exp 



/3 2 ^A 



X(£) 



(58) 



and therefore, at i = £*, where X(£*) = 2n, we have e e ~ LqX^ 3 ^ exp[27r/(3Ao)]. Conse- 
quently, C(t, A) w (7(0, A) w C A W , with w = /5 3 //3 2 2 . 

Within the approximations leading to Eqs. (|53|) and (|5^) , we note that the expressions 
in the exponentials and their first derivatives are continuous at Xq = to, up to terms of order 
0(Ao/27r). In comparing our results with the experiment, we shall use these asymptotic 
expressions all the way to the line A = t . 
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V. THE PHASE BOUNDARY FOR WEAKLY COUPLED PLANES 



The 3D transition temperature T^(x), as function of the defect concentration x, of a 
system consisting of weakly coupled planes may be deduced from the relation 

ag D (t N ,\)~l. (59) 

[Note that A is proportional to x, cf Eq. (0).] The parameter a can be generated by an 
interplane exchange, J±/J ~ J±/2irp s , or some in-plane spin anisotropy. As stated in the 
Introduction, this procedure gives excellent estimates for T/y. || 

To obtain the critical line T^{x) we proceed as follows. Using fl52D in the relation (|59|) 
we find 

l-^ = ^lnKL C(t,A)) 2 ]- 1 . (60) 

At low temperature, i.e., for A > to, we have X(£*) = 2ti and therefore Eq. ( |60"D is almost 
independent of t. It thus yields a critical value for the initial value of the variance, A c , 
and hence a critical concentration, x c , above which there is no antiferromagnetic long-range 
order at any temperature 

2tt 

^-l + flnKLo^CA,)) 2 ]- 1 ' (61) 
with A c = Ax c . In fact, £ 2 d is expected to be practically independent of temperature (except 
for the very weak dependence of the prefactor C) for a range of values, A > to, as given 
in Eq. (|53|). Therefore, the critical line T^(x) is expected to be practically vertical for 
tj^(x) < Ao = Ax. Below this line, one might expect some range in which spin glass and 
antiferromagnetism co-exist, down to a Gabay-Toulouse like line. |17j To obtain this region 
one would need to also consider the 3D boundaries of the spin glass phase, and this is beyond 
the scope of the present paper. In any case, the AFM ordering will persist up to the line 

At smaller defect concentrations, or at higher temperatures, i.e., for 2-n ^> t > A , Eqs. 
(0) and © give 



Ap _ f _ 3A -i 1/3 

t N (x) ~ <■ ~ t N (0)i ' 



= ^InlaCLoCMOj.O)) 2 ]- 1 , (62) 
where ijv(O) is the Neel temperature of the pure antiferromagnet. The lowest order of this 



expression agrees with the results of Ref. [14], obtained by an annealed average. 



VI. ANNEALED AVERAGING 

As stated above, the cuprates require some mixed annealed-quenched averaging. We 
start be reviewing a simple version of GI's analysis. Jl4| In their approach, and a(r^) are 



treated as quenched variables, with averages given by Eq. (||), while m(r^) is treated as 
annealed. Thus, Eq. @ is replaced by 

[/v(ri)/j>(r 2 )] = 5 ij 5(ri 2 )A AtI/ (r 1 ), (63) 

with 

A^(r) = mil {r)m u {r)M 2 x/d. (64) 

Initially, Ti contains no interactions among the dipole moments {m(r)}. However, the 
RG iterations generate the dipole-dipole interaction, as given by Eq. (|28|) . This interaction 
is mediated via the canted background AFM moments. Treating this interaction as a small 
perturbation, to lowest order, we can next integrate the dipole moments out of the partition 
function, using the annealed averaging (m M (ri)m„(r 2 )) = ^5(ri 2 )/A/', so that 

(A^> = 5^ A = 5^M 2 x/(dAf). (65) 

Note that to this leading order, A = A! GI wrote down a more general form for this thermal 
average, involving the susceptibility which results from the quadratic coupling in Tidd- This 
reduces to the above expression at lowest order. 

Up to Eq. (P^D, we have performed no averaging. In the annealed case, the recursion 
relation are derived from the same equations, using the averages as listed above. The 
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resulting recursion relation for 1/t turns out to be the same as Eq. fl34|), with A replaced 
by A. In contrast, the averaging over m gives no contributions to the renormalization of fj, 
since all the generated terms which are linear in B- 1 involve odd powers of the u^s and the 
gf s, and therefore odd powers of the m^'s. All of these vanish upon the annealed averaging 
over the m^'s. Thus, we end up with 

(£)' = *"(£)■ («) 

At 2D, Eq. (^) implies that A is unrenormalized. The solution of the recursion relation 
for t then yields 

* ln l + 2AA(0) 

A 1 + 2K(t{€) ' y ' 

Assuming that A <C 2tt, and integrating up to t(t) = 2tt, yields 

^^(l-A), (68) 

which agrees to leading order with the quenched result (|5^). 

This annealed averaging is legitimate as long as the renormalized distance beteen the 
impurities remains large, so that the dipole-dipole interaction which is generated during the 
iterations remains small. Note that the range of Tidd is 6 = e . Thus, if e is small compared 
to the inter-impurity distance x~ x l A : , then we can still treat the dipoles as independent degrees 
of freedom, ignore the interaction between them and continue the above annealed calculation. 
However, if e e * becomes larger than (but of the order of) x~ l l d then each renormalized cell 
contains more than one impurity, and the interaction between them comes into play. Since 
this interaction decays as 1/r , the dipoles behave like a spin glass which is at its lower 
critical dimension. |T7[] Since the dipole-dipole energy is of order E d d = Ap s , we expect 
the dipole moments to develop spin-glassy correlations, with a correlation length £ sg which 
grows exponentially in Edd/T. Since we showed that £20 remains finite at all T, we expect 
that for T E d d one has £ sg 3> the moments inside a renormalized cell (at distances 
smaller than ^d) freeze randomly, with the effective Edwards-Anderson order parameter 
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Q, and we can switch to our quenched calculation. We thus choose to perform a prefacing 
annealed renormalization, up to e e ° = L ~ x~ x l d . For larger t we assume that the dipole 
moments are frozen at low T, and we switch to the quenched analysis of Sees. IV and V. 

VII. COMPARISON WITH EXPERIMENTS ON DOPED CUPRATES 

Given the discussion in the previous section, we adopt the following strategy: We start 
by fitting A = X/x from the t-dependence of at high T and small x. As stated, this 
dependence is the same for both types of averages. Given A, we next fit the x dependence of 
^2D in the limit of very low T, when is T- independent. This behavior should certainly be 
described by our quenched theory. The fit determines the pre-exponential factors. Finally, 
we use the results (without any further adjustments) to calculate the phase diagram, T N (x). 

We begin with the temperature dependence of the correlation length at low concentra- 
tions, x < x c . Data taken on a sample of La 2 Cu04 +< 5, with T N = 90K, j|] show a practically 
linear dependence of (i/27r) ln(£/C) on 1/t, in agreement with both our quenched result fl54|) 
and our annealed result (|68|). For the coefficient C in this fit we used C = 1.92A, derived 
from Eq. © with p s = 24meV (See Ref. |]) and L = 1. The slope is A = A = 0.29(1) 
(see Fig. 1, and also Ref. [|15|]). To estimate the value of x for this sample, we follow Refs. 
Q and |29[] and approximate the line T^(x) by the straight line T^(x) ~ 325 — 16250a;, 
which extrapolates to x = 0.02 as Tjy is extrapolated to zero. Using this approximation, we 
find that Tjv = 90K at x = 0.0145. Thus, A = X/x ~ 20. Although we have some problems 
with this linear extrapolation (see below), the value of x cannot be larger than x c ~ 0.02, 
so the uncertainty in A is not more than 30%. Furthermore, although A might have a weak 
dependence on T and on x, this could most probably be absorbed in this error estimate. We 
thus use this estimate A = 20 in everything that follows. 

Keimer et al. measured the temperature dependence of the correlation length for three 
magnetically disordered samples of L^-xSr^CuO^ with x ~ 0.02, 0.03, and 0.04. The error 
in x is less than ~ 0.005. It is believed [ffl that the hole concentration is about that of the 
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Sr ions. It was found that at low temperatures £ does not depend on T, and falls with the 
increase of x faster than a power law. The £(T = 0) data is depicted in Fig. 2, together with 
the value of £ cited in Ref. [12|] for Lai.gsBao.osCuO^ The figure also shows the theoretical 
values of £, calculated from Eq. fl53|), with L C(0, Ao) = Cq\q , Ao = 20a;, Co = 2.8A and 
a/ = 0.8. (In 2D, we now have u' = uo — 1/2, from the x dependence of L ). We also 
reproduce in Fig. 2 the numerical results obtained in Ref. [^7]]. These authors computed 
the effect of holes localized on the oxygen atoms in the CuC>2 plane. Their results, marked 
by '+' in Fig. 2, are in agreement with our calculation. 

We mention in passing that the correlation length at low temperature was also measured 



for x > 0.05. [|30| These data were not included in our comparison, mainly because for these 
concentrations the holes are probably mobile. 

We next turn to the temperature dependence of £ at x > x c . Figure 3 exhibits the 
calculated temperature dependence of £, for several concentrations x. The curves were 



found from Eq. fl5"4|) with C = L C = 1.92A (Ao < t , as in Fig. 1 discussed above), and 
from Eq. Q53[), with C{x) = 2.8X° 8 A (A > t ). We have also used 2vrp s = 150 meV. The 
theoretical lines are for x= 0.0225, 0.029 and 0.036 (instead of the experimental values 0.02, 
0.03, and 0.04 given by Keimer et al.). The values chosen are within the experimental error. 
H] Since the prefactors used in our fits differ in the limits A <C t Q and A t Q , the high- 
and low-temperature portions of the plots do not match at Ao = to- There the two segments 
are connected by dotted lines. The dashed lines in Fig. 3 represent the heuristic expression 

C 1 (T,x) = C 1 (0,x)+C 1 {T,0), (69) 

used in Ref. [Q to fit their data. It seems that the heuristic expression works as well; however, 
it seems that it has no theoretical basis. 

The theoretical curves in Fig. 3 agree with the experimental results at temperatures 
lower than ~ 350 — 400K. At higher temperatures the calculated values of £ are smaller 
than the measured ones. Perhaps, at such high temperatures thermal fluctuations come into 
play, causing a decrease of Q from 1/jV to lower values. It is also possible that at high T 
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the holes are more mobile, so that one should average over more than one bond per hole, 
thus reducing the effective dipole moment. At temperatures lower than 200 — 250 K the 
correlation length does not depend on the temperature up to exponentially small terms, of 
the order of £(0, T) -1 . [|J This nontrivial property of the correlation length is reproduced 
well by our calculation. 

Given the above values for Co and u, and the value a = 10~ 4 from Ref. ||, we solved 
Eq. (|6l| ) and found the critical value of A to be A c = 0.366. With A = 20, the critical 
concentration, x c , for the disappearance of the long-range order at T = is found to be 
x c = 0.0183, in very good agreement with the experimental value x c = 0.0175 from Ref. || 
(but in disagreement with Ref. ||, which gives x c = .027). 

Using the above parameters, i.e. A = 20, a = 10 -4 , and approximating the prefactor 
by a constant, C(t,x) « C(0,x c ) = (7(0,0.0183) = 1.26A, Fig. 4 shows the theoretical 
concentration dependence of T/v, calculated from Eqs. ( |62^ ) for t > X Q and flBTD for \ > t , 
with no further fitting of the parameters. (The line t = A is also depicted in the figure.) At 
high temperatures the prefactor is larger, 1.9A. However, the effect of this difference on T/v 
is small, since T/v depends on the prefactor only logarithmically. At small concentrations, 
our theoretical T^{x) decreases with the increase of x linearly with the rate 55 K/% (based 
on the value of A as determined from the £ data). At x — x c , T/v(x c ) is roughly equal to 
7jv(0)/3. Then T/v abruptly falls to zero, as our calculation finds that the correlation length 
is independent of the temperature for t < A . 

Figure 4 includes the results of several experiments. It is seen that those of Ref. |29| 
seem to disagree with our phase diagram: the data fall linearly with a slope of 162 K/%, 
(larger by about a factor of 3 than our theoretical value, which was extracted from the data 
for £), extrapolating to x — .02 without the jump in T^. However, we should note that Chen 
et al. |29j determined x for their O-doped samples from Hall effect data. In Ref. T/v 
and the Hall density of holes were measured for a Sr doped sample. The Sr concentration 
in the melt was 0.0022, while the Hall measurements gave a smaller hole density, 0.0016. 
Figure 4 shows that in this case the experimental point (with x = 0.0022) is closer to the 
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theoretical curve than the data from Ref. ||29|| . Saylor and Hohenemser || measured Tjy(x) 



in Sr-doped samples of lanthanum cuprate. Although their Tjv(O) = 317 K was somewhat 
less than in the best samples of Ref. |4]], their Tn(x) decreased linearly with the rate 90 K/% 
till x = 0.015. In the region between x =0.015 and 0.018, T N fell from « 180 K to 20 K. 
This behavior of T^(x) is close to our phase diagram, Fig. 4. 

Very recently, Wakimoto et al. jT(| measured the phase diagram of oxygen doped 
Lai.9 5 Bio.o5Cu04, and their Tn(x) agrees qualitatively with our theory: it falls almost lin- 
early down to T/v (0.012) m 160K, and then drops sharply towards x c ~ 0.015. Both this 
small value of x c and the low-T value of the correlation length near x c , as measured in Ref. 
T0|1 , are consistent with our calculations, with A m 30. 



VIII. DISCUSSION 

Most of this paper was devoted to a detailed description of the theory for the effects of 
FM bonds on AFM correlations in doped antiferromagnets. We hope that this will stimulate 
some material research onto such systems. It is also hoped that this exact solutions could 
be used to study more systems with correlated random fields. 

In the previous section we discussed some fits of our theory (with both annealed and 
quenched averaging) to data from doped lamellar oxides. As stated in the Introduction, these 
fits are based on various assumptions, and thus deserve some further discussion. Hopefully, 
the reader is already convinced that our fits give a nice interpolation between the annealed 
and quenched calculations, for high and low temperatures respectively. We now discuss 
possible problems in using our model of localized holes for these systems. 

First of all, we note that experiments on both O-doped and Sr-doped La2Cu04 indicate 
a localization length of the hole of order of two lattice constants. p9| , |3Tf Furthermore, for 
Sr-doping, the hole moves around the center of a Cu-plaquette, and this might imply a 



quadrupolar behavior, rather than a dipolar one, for the spin canting. [15| However, there 



are arguments in the literature, claiming that holes which move around such a center 
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create Skyrmion-type distortions of the AFM background, decaying as 1/r. Thus, our 
calculation may remain qualitatively correct, although we might have to use a different 
expression for our parameter A. We might also add that to linear order in x one effectively 
considers only single impurities, and therefore mobile and static holes will have the same 
(classical) effects. Since A is fitted from the experiments anyway, it could result from any 
of the above cases. 

Another alternative is to abandon the localized picture altogether. The effects of mobile 
holes on the AFM correlations were considered in many papers. One approach |$3| concluded 
that mobile holes drive the system to the quantum disordered state. Indeed, we solved the 
classical NLaM, and thus ignored quantum fluctuations. However, for the undoped systems 
the classical expressions described the exponential factor of the correlation length rather 
well. [|T9|,[4|| We have checked that integration of the quantum fluctuations, similar to that 
done in Ref. |]I9|| , only renormalizes the initial parameters of the effective classical model 
also in the random case. 

Another approach proposed that the AFM order is destroyed by the segregation of 
holes into striped phases. In view of the experiments metioned above, which show that the 
holes are localized (and the conductivity obeys the variable range hopping law) at low T 
and low x, we believe that these alternative approaches should be used only for x > 5%, or 
at higher T. Indeed, to our knowledge there is little experimental evidence that stripes exist 
at x < 5%, or that quantum fluctuations were observed at low T in that range. However, 
we cannot rule out some crossover from our localized to the mobile behavior at the metal- 
insulator transition (around x = 5% and at higher T). In any case, it would be nice to 
find more experimental tests that would distinguish between the various scenarios for the 
cuprates. 
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APPENDIX A: GAUGE TRANSFORMATION OF THE FIELDS <p^ 

In order to eliminate the functions Af u (r) from the expression for the derivatives c^n 
[Eq. (|I5D 1 we introduce the transformation 



0M = E r w*An ( A1 ) 

w 

where the coefficients T Wl are determined by 

djTu Ul = >^ Af M2 T M2W1 , (A2) 

with 

T ulM = — (A3) 

Note that the sums run from 1 to M — 1. Assuming that Af 11 is independent of x,, the 
solution of Eqs. ( |A2| ) and ( |A3| ) reads T MW = T° exp(— iKXi), where k is a real eigenvalue 



of the Hermitian matrix iA{. Deviations from this approximation naturally involve higher 
derivatives of A^ 1 which are related to higher powers of the gradients in TC. These are 
strongly irrelevant. 

Inserting the above relations into Eq. (pl|) we obtain 



+ E [B?ejl - 2 + ditp E T V(1 e v ] . (A4) 

Therefore, defining 
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e M = $>T^, B* = Y,BiT»i» (A5) 



Eq. ( |A4j) takes the form of Eq. (p~5|) without the terms, with M , and replaced 



by M , e M and £>f , respectively. For brevity, we have omitted in the subsequent derivations 
the superscript ~ It is seen that the gauge transformation can be regarded as a rotation of 
the base vectors e M , and reflects the arbitrariness of their choice. p4| 

APPENDIX B: THE PHASE BOUNDARY AT 2+e DIMENSIONS 

For completeness, we present here the solution of the recursion relations (|46|) in d = 2 + e 
dimensions, and obtain the critical line t^(x), where x is the defect concentration. 

The recursion relations ( [46|) have four fixed points in the [A, t] plane: [0, 0], [0, 27re/(A/" — 
2)], [27re/7V,0], and [2vre/(27V - 3), 2vre/(27V - 3)]. For M > 3/2, the first point is stable, 
the last one is doubly unstable and the other two are unstable in one direction and stable 
in the other. 

The solution of Eqs. ( [46]) for i > £ and Heisenberg spins (A/ - = 3) reads 

L 2ire j 

l-X(£)/t(£) r A(£)-2 7 re/3 1 i/3_ 



1 - Ao/to L A - 27re/3 

The phase boundary is identified as the line which separates the flow to the origin (in 
the ordered phase) from the flow to infinity (in the disordered phase). It is easy to check 
that for to > Ao, the initial point [Ao, t ] will flow to the "pure" fixed point [0, 2ire] when 

A /, 3A xi/3 



27T6 

< Ao < A c = — , (B2) 



where the second of Eqs. ( pi|) has been used. On the other hand, using again that equation, 
we find that for t < X the solution flows to the "random" fixed point, [27re/3, 0], provided 
that 
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2ne 

\ = \ c = A c >t >0. (B3) 



The two portions of the critical line, Eqs. ( |B2| ) and ( |B3| ), are separated by the multicritical 



point at [27re/3, 27re/3]. It is interesting to note that also here, like the behavior found for 
t]y(x) for weakly coupled planes, there is a vertical section of the phase boundary. 

We thus conclude that for to > Ao the randomness is irrelevant, and the "pure" critical 
behavior (which has a correlation length exponent v = 1/e) remains stable. However, 
finite values of A yield a correction to this behavior, with exponent — e (or, more generally, 
— e/ (A/"— 2)). The correlation length is obtained from Eq. (0), using the matching conditions 
(H). One finds 

A — A c 

For t > Ao the iterations are stopped at t(£*) ~ 2ti. In the range X(£*) < A c we can find 
\{t) by considering the difference \o/t N — \o/t , using Eqs. ( [BT| ) and flB2"|). We find 

tf-^-V^d-^)" 3 , (B5) 

from which it follows that 

£~e e ' ~(t -Mx))- 1/£ x 

[1 + Sar(t - ^(ar)) + ...], (B6) 

where i? is a constant and the term associated with it comes from corrections of order Aoe _r . 
For A > to, A(r) = 2ir, so that 

£~ e r - (A 1 - A; 1 )- 1 ^. (B7) 
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FIGURES 
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FIG. 2. Dependence of at T = on x. The empty 4 and full 12 circles indicate experiments, 
+'s show numerical simulation 27 data. The solid line represents 2.8A°' 8 Aexp(27r/3A), with A = 20x. 
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FIG. 3. Dependence of £21? (i, A) on T for several concentrations. Symbols are from experiments: 
empty circles for x = 0.02, full circles for x = 0.03, +'s for x = 0.04, all from Ref. 4. Full lines 
show results from Fig. 2 (for low T) or Eq. (54), with C = L C = 1.921 (for high T). Dotted 
lines interpolate between these low- and high-T theories. Dashed lines correspond to Eq. (69). 
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Fig. 4 
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FIG. 4. Tjv(x)/T/v(0) versus x. Full line is theory, and the points are from experiments: full 
circles from Ref. 6 , +'s from Ref. 7, empty circles from Refs. 29 and 31. Dashed line indicates 
Aq = to- 
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